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We present a preliminary measurement of the Standard Model CP violation parameter sin 2</>i using 
the Belle detector. A 6.2 fb" 1 sample of events produced by the KEKB asymmetric e + e~ collider 
operating at the T(4S) resonance is used. One neutral B meson is fully reconstructed via its decay to 
a CP eigenstate: J/ipK$, ip{2S)K$, Xcl^Si J /M>K-L or J/ip-K . The flavor of the accompanying B is 
identified mainly from the charges of high momentum lcptons or kaons among its decay products. The 
time interval between the two decays is determined from the distance between the decay vertices. A 
maximum likelihood fitting method is used to extract sin 20i from the asymmetry in the time interval 
distribution. We report a preliminary result of 

sin 2<^i = 0.45±°;«(stat)±°;^(syst). 



1 Introduction 

The Belle experiment at the KEKB asym- 
metric energy e + e~ B-meson factory recently 
completed a successful first year of operation. 
The KEKB luminosity which was about 
10 31 cm~ 2 s at the time of startup in June 
1999, reached a level of 2 x 10 33 cm" 2 s _1 by 
July 2000. The Belle detector collected a to- 
tal integrated luminosity of 6.8 fb _1 , 6.2 fb _1 
on the T(45) resonance and 0.6 fb _1 off res- 
onance. We use this data sample to make 
a preliminary measurement of the Standard 
Model CP violation parameter sin20i using 
B° d J/^K s ,^(2S)Ks,XaiK s ,J/^7r° and 
J/4>Kl decays^ 

The Standard Model predicts a, CP vio- 
lation through a mechanism of Spontaneous 
Symmetry Breaking of the electroweak sym- 
metry that results in the Cabibbo-Kobayashi- 
Maskawa (CKM) quark mixing matrix. In 
systems of mesons containing a b quark, CP 

"Throughout this paper, when a mode is quoted the 
inclusion of a charge conjugate mode is implied unless 
otherwise stated. 



violating effects are expected to be large. 
The interference between the direct B® — > 
fcp decay amplitude and the mixing- induced 
B® — > Bd —> fcp decay amplitude, where 
fcp is a CP eigenstate to which both B® and 
Bd can decay, gives rise to an asymmetry in 
the time-dependent decay rate: 

A u\ = dN/dt(Wl a ^fcp)-dN/dt(B a t=0 ^f C p) 
^ ' dN/dt(B" =0 ->fcp)+dN/dt(B° t=0 ->fcp) 

= —rjf sin 20i sin Am^i, 

Jl) 

where t is the proper time, dN/dt(B t=0 
(5Lo) fcp) is the decay rate for a B°(B°) 
produced at t — to decay to fcp at time 
t, rjf is a CP-eigenvalue of fcp, (Vf = ~ 1 
for J/ipK s ,ip(2S)K s and XaK s , and +1 
for J/ipir° and J/ipKr,), Am<j is the mass 
difference between two B° mass eigenstatea. 
and 0i is one of the three internal anglesEl 
of the CKM Unitarity Triangle, defined as 
^Tr-arg^g^). 

In T(4S) decays, B° and 23° mesons are 
pair-produced and remain in a coherent p- 
statc until one of them decays. The decay 



of one B° meson to a final state /i at time 
t\ projects the accompanying B meson onto 
an orthogonal state at that time; this meson 
then propagates in time and decays to f% at 
time t2- CP violations can be measured if 
one of B mesons decays to a tagging state, 
ftag, i-e. a final state unique to B° or B , at 
time ttag and the other decays to an fcp state 
at time tcp- A time-dependent asymmetry, 
A(At), which is obtained by replacing time 
t in (|l|) with the proper time interval At = 
tcp ~ ttag, can then be observed in T(45) 
decays. Because the B°B° pair is produced 
nearly at rest in the T(4S) center of mass 
system (cms) , At can be determined from the 
distance between fc p and f tag decay vertices 
in the boost (z) direction, zqp and z tag , and 
At ~ Az/(3jc, where ^7 is a Lorentz boost 
factor of the T(45) restframe (equal to 0.425 
at KEKB). 

This asymmetry is diluted by experimen- 
tal factors including background to the re- 
constructed fcp states, the fraction of events 
where the flavor of the B meson is incorrectly 
tagged (lu), and the resolution of the decay 
vertex determination (d res ): 

-^-observed = I J _|_ BIS ~~ 2w)ct res }^4 = DA. 

(2) 

where B/S is the ratio of background to sig- 
nal and D(< 1) is the "dilution factor." The 
statistical error of sin 2(/>i is inversely propor- 
tional to D: (5 sin 20i = -^==j^. 

2 The Belle Detector 

Belle is a large-solid-angle magnetic spec- 
trometer 0. Charged particle tracking is pro- 
vided by a three-layer, double-sided silicon 
vertex detector (SVD) and a small-cell cylin- 
drical drift chamber (CDC) consisting of 50 
layers of anode wires, 18 of which are inclined 
at small angles. The tracking system is situ- 
ated in a 1.5 T solenoidal field. The charged 

^Background in (feh is assumed to have no asymmetry. 



particle acceptance is 17° < 6 < 150°, where 
9 is the polar angle in the laboratory frame 
with respect to the beam axis; the corre- 
sponding cms acceptance is ~ 92% of the 
full solid angle. The impact parameter res- 
olutions are a 2 r<j) = (21) 2 + ( pf3 ™ 3/20 f P-m 
in the plane perpendicular to the beam axis, 
and a 2 = (39) 2 + ( pf3s ^ /2 e ) 2 Aim along the 
beam direction, where p is the momentum 
measured in GeV/c and (3 is the velocity di- 
vided by c. The transverse momentum reso- 
lution is (cr pt / P t) 2 = (0.0019p t ) 2 + (0.0034) 2 . 
Charged hadron identification is provided 
by dE/dx measurements in the CDC, aero- 
gel Cherenkov counters (ACC) and a bar- 
rel of 128 time-of- flight scintillation coun- 
ters (TOF). The dE/dx measurements have 
a resolution for hadron tracks of a{dE/dx) = 
6.9% and are useful for tt/K separation for 
p < 0.8 GeV/c. The TOF system has a 
time resolution of 95 ps (rms) and provides 
tt/K separation for p < 1.5 GeV/c. The 
indices of refraction of the ACC elements 
vary with polar angle from 1.01 to 1.03 to 
match the kinematics of the asymmetric en- 
ergy collisions and provide tt/K separation 
for 1.5 GeV/c < p < 3.5 GeV/c. Parti- 
cle identification probabilities are calculated 
from the combined response of the three sys- 
tems. The efficiency for K ± is ~ 80% with a 
charged pion fake rate of ~ 10% for all mo- 
menta up to 3.5 GeV/c. 

An array of 8736 CsI(Tl) crystals pro- 
vides electromagnetic calorimetry that cov- 
ers the same solid angle as the charged 
particle tracking system. The photon en- 
ergy resolution, estimated from beam tests, 
is {cte/E) 2 = (0.013) 2 + (O.OOO7/£0 2 + 
(O.OO8/.E 1 / 4 ) 2 , where E is measured in GeV. 
Neutral pions are detected via their decay 
to 77. The 7T° mass resolution varies slowly 
with energy, averaging a m — 4.9 MeV/c 2 . 
With a ±3cr mass selection requirement, the 
overall detection efficiency, including geomet- 
ric acceptance, for ir°s from BB events is 
40%. Electron identification is based on a 



combination of the CDC dE/dx information, 
the response of the ACC, and the position, 
shape and energy deposit of its associated Csl 
shower. The electron identification efficiency 
is above 90% for p > 1.0 GeV/c with a pion 
fake rate that is below 0.5%. 

The magnetic field is returned via an 
instrumented iron yoke consisting of alter- 
nating layers of resistive plate counters and 
4.7 cm thick iron plates. The total iron thick- 
ness of 65.8 cm plus the material of the Csl 
calorimeter corresponds to 4.7 nuclear inter- 
action lengths at normal incidence. This sys- 
tem, called the KLM, detects muons and Kl 
mesons in the region of 20° < 8 < 155°. 
The overall muon identification efficiency is 
above 90% for p > 1 GeV/c tracks detected 
in the CDC; a pion fake rate is below 2%. 
Kl mesons are identified by the presence of 
the KLM hits originating from hadronic in- 
teractions of the Kl in the Csl and/or iron. 
If there are Csl hits associated with a can- 
didate Kl, its direction is determined from 
the energy- weighted center of gravity of the 
Csl hits alone. Otherwise, the direction is 
determined from the average position of the 
associated KLM hits. The angular resolution 
of the Kl direction is estimated to be ~ 1.5° 
and ~ 3° with and without associated Csl 
hits. 

3 Selection of B° Decays to CP 
Eigenstates 

We reconstruct B° decays to the following 
CP cigcnstateJ: J/ipKs, ^{2S)K S , XciK s 
for CP — 1 states and J/ipn , J/i/jKl for 
CP+1 states. The B°B° hadronic events are 
selected by requiring (i) at least three tracks 
with a minimum p t of 0.1 GeV/c originating 
within 2.0 cm and 4.0 cm of the run-by-run 
average interaction point (IP) in the plane 
transverse to the beam axis {xy plane) and 
along the beam axis (z axis), respectively; 
(ii) at least two neutral energy clusters with 
a minimum energy of 0.1 GeV in the barrel 
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Figure 1. The invariant mass distributions for (a) 
J/4> -> H+H~, (b) J/ip -> e+e~ and (c) i/>(2S) -> 
fjb + y,~ , e+e - . 



Csl calorimeter; (iii) a sum of all Csl clus- 
ter energies that is between 10 % and 80 % 
of the cms energy (B cms ); (iv) the total visi- 
ble (charged and neutral) energy greater than 
20 % of £ cms ; (v) \ J2pT s \ les s than 50 % of 
E cms /c, where pl ms is the z component of the 
momentum calculated in the cms frame; (vi) 
a reconstructed event vertex within 1.5 cm 
and 3.5 cm of the IP in the xy plane and 
along the z axis, respectively. In addition, we 
apply an event topology cut, H 2 /H < 0.5, 
where H% and Hq are the second and zeroth 
Fox- Wolfram moments, to reject continuum 
background. 

The J/tp and ip(2S) mesons are recon- 
structed via their decays to and e + e~. 



Dimuon candidates are oppositely charged 
track pairs where at least one track is pos- 
itively identified as a muon by the KLM sys- 
tem and the other is either positively iden- 
tified as a muon or has a Csl energy de- 
posit that is consistent with that of a min- 
imum ionizing particle. Similarly, dielectron 
candidates are oppositely charged track pairs 
where at least one track is well identified 
as an electron and the other track satisfies 
at least the dE/dx or the Csl E/p electron 
identification requirements. For dielectron 
candidates, we correct for final state radia- 
tion or bremsstrahlung in the inner parts of 
the detector by including the four-momentum 
of every photon detected within 0.05 ra- 
dian of the original electron direction in the 
e + e~ invariant mass calculation. The in- 
variant mass distributions for J/ip — > 
J/tp — ► e + e~ and tp(2S) — ► /i + /i~,e + e~ 
are shown in Figs. 0(a), (b) and (c), respec- 
tively. The ip(2S) is also reconstructed via its 
J/ipTr + TT~ decay and the Xci is reconstructed 
through its decay to J/ipj. Figures ^ (a) and 
(b) show the mass difference distributions of 
Mi+i-K+x- ~ M e+e - and M i+e -^ - M e+t -. 

Candidate K$ — > n + ir~ decays are op- 
positely charged track pairs that have an in- 
variant mass between 482 and 514 MeV/c 2 , 
which corresponds to the ±3cr around the K$ 
mass peak. The Ks — > 7r 7r° decay mode 
is also used for the J/ipKs channel. These 
are selected among photons with a minimum 
energy of 50 MeV and 200 MeV in the bar- 
rel and endcap regions, respectively, by re- 
quiring, assuming Ks decayed at the IP, (i) 
a minimum ir° momentum of 100 MeV/c; 
(ii) 118 < M 77 < 150 MeV/c 2 ; and (iii) 
300 < M^o < 1000 MeV/c 2 . For each 
candidate, we determine the most probable 
K$ decay point by minimizing the sum of 
the \ 2 values from constraining each photon 
pair to 7T° invariant mass while varying the 
Ks decay point along the Ks flight direction 
defined by the sum of four photon momenta 
and the IP. We then recalculate the invariant 
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Figure 2. the mass difference distributions of (a) 
Mi+i-v+n- - M l+l - and (b) M e+l - 1 - M l+l - . 



masses of the photon pairs and the 7r°s and 
require the recalculated Ks mass to be be- 
tween 470 and 520 MeV/c 2 . For the J/tpir 
mode, the tt° candidates are selected from 
photons with a minimum energy of 100 MeV. 

To identify reconstructed B meson de- 
cays we use the beam-constrained mass 
M beam = \/ £f e am _ Pb and tne energy 
difference AE = Eb — E beam , where 
E beam is E cms /2 and p B and E B are 



the B candidate three-momentum and en- 
ergy calculated in the cms. Figure || 
shows the M beam distribution of the com- 
bined B -> J/ipK s (TT + Tr-), J/^K s (n n ), 
t/;(2S)Ks(n+Tr-),XciKs(Tr + Tr~), and J/^tt 
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Figure 3. The beam-constrained mass distribution 
of sum of B -> J/i/.E's(i" + 'r"), J /iI>Ks{it -k ), 
4>{2S)Kg(ir~l~ir~),XciKg(iT + iT~), and J/ipw events. 



samples, after the imposition of a 3.5cr cut 
on | AS | (±40 MeV for modes with K s -> 
7r + 7r _ and ±100 MeV for modes containing 
7r°s). The -B meson signal region is defined 
as | Mb 

earn < a in 

> | < 0.01 GeV/c 2 , 
where < Mb eam > is the mean value of ob- 
served Mbeam- Table lists the number of 
signal candidates (N ev ) and the backgrounds 
(Nbkgd) determined from extrapolating the 
event rate in the non-signal AE vs Mbeam 
region into the signal region, and from the 
full Monte Carlo (MC) simulation results. 

The B° — > J/tpKL event candidates 
are selected by requiring the J/ip momen- 
tum and the Kl direction to be consistent 
with two-body-decay kinematics. After re- 



Table 1. Summary of the reconstructed CP eigen- 
states 




Mode 


N ev 


Nbkgd 


J /%jjKs{'K + 'K~) 


70 


3.4 


J/^K s (tt°tt ) 


4 


0.3 


ip(2S)K s (Tr+n-) 


5 


0.2 


ip{2S)(J/ifjir+TT-)K s (Tr+n-) 


8 


0.6 


Xci(iJ/^)K s (tt + tt^) 


5 


0.75 




10 


1 


Total 


102 


6.25 



» 1.0 

P B (GeV/c) 

Figure 4. The p* B distribution with the results of the 
fit. The upper solid line is a sum of the signal and 
background. The total background (lower solid line) 
is divided into those coming from J/ipK*° (K^tt ) 
and Ji\) + non-resonant K^tt (above the dotted line) 
and those coming from all other sources (below the 
dotted line). 



quiring the J/ip cms momentum to be be- 
tween 1.42 and 2.0 GeV/c, we calculate the 
cms momentum of the B meson, p* Bl which 
should be equal to ~ 0.34 GeV/c for a true 
event. Figure ^ shows the p* B distribution. 
Also shown are expected distributions of sig- 
nal and background derived from the full MC 
simulation studies. The background is found 
to be dominated by B — > J/tpX events in- 
cluding B -> J/ipK*°(K L n°) and B -> J/ip 
+ non-resonant K^tt , which are a mixture 
of CP + 1 and CP — 1 eigenstates. There are 
102 J/tpKi, candidates in the signal region 
defined as 0.2 < p% < 0.45 GeV/c. By fitting 
the data with the expected shapes, we find 
48 background events in the signal region, of 
which 8 events were from J/iPK* (Kltt ) + 
J/ip non-resonant Kjjn . 

4 Flavor Tagging 

Each event with a CP eigenstate is exam- 
ined to see if the rest of the event, defined as 
the tag side {ftag), contains a signature spe- 
cific to B° or B . Our tagging methods are 



based on the correlation between the flavor of 
the decaying B mesons and the charge of a 
prompt leptons in b — > cvi decays, the charge 
of Kaon originating from b — > c — > s decays, 
or the charge of ir from i? — > D*(—^ irD)£v 
decays. A B°(B ) flavor for f tag indicates 
that fcp was in B (B°) state at At = 0. 
We applied the following four tagging meth- 
ods in descending order: if the event failed 
the method (1), we tested the event with the 
method (2), etc. 

1. High momentum lepton: If f tag contains 
a lepton = e ± or /i^) with cms 
momentum p* > 1.1 GeV/c, we assign 
f tag = B°(B°) for £+(£-). 

2. Charged Kaon: If f tag contains no high 
momentum £^ , the sum of the charges 
of all identified Kaons, Qk, in ftag is 
determined. We assign f tag = B°(B ) if 
Qk > (Qk < 0). If Q K = 0, the event 
fails this method. 

3. Medium momentum lepton: If f tag con- 
tains an identified lepton in the cms mo- 
mentum range 0.6 <p}< 1.1 GeV/c, we 
use the cms missing momentum [p* miss ) 
as an approximation of the v cms mo- 
mentum. If p*f +p* miss > 2.0 GeV/c, we 
assume f tag is from b — > ci^ decay and 
assign its flavor based on the charge of £ 
as in the method (1). 

4. Soft pion: If f tag contains a low mo- 
mentum (p* < 200 MeV/c) charged 
track consistent with being a 7r from 
the D* — > D7r decay chain, we assign 

/ tog = for TT" (7T+). 

The efficiency (e) and the wrong tag fraction 
(u>) are determined using a sample of exclu- 
sively reconstructed B a — > D*~ (D~)£ + v de- 
cays which are self-tagging decay modes, and 
the full MC simulation. Because of B° — B 
mixing, the probabilities of finding the oppo- 
site flavor (OF) and same flavor (SF) neutral 
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Figure 5. The B° — B mixing amplitude as a func- 
tion of the proper-time interval At of two neutral 
B mesons obtained using (a) B° —* D*£u decays 
and (b) B° — » Dlv decays. Also shown are the re- 
sults of the fit (solid lines). We obtained Amj = 
0.488 ±0.026 pb _1 . 



Table 2. Tagging efficiency (e) and wrong tag fraction 



Method 


e 


uo 


High p* £ 
K± 

Med. p* £ 

Soft TT 


0.0142 ±0.021 
0.279 ±0.042 
0.029 ±0.015 
0.070 ±0.035 


0.071 ±0.045 
0.199 ±0.070 
0.29 ±0.15 
0.34 ±0.15 



B meson pairs are 

PoF(At) oc 1 + (1 - 2w) cos(Am d At) 
P SF (At) oc 1 - (1 - 2w) cos(Am d Ai). [ ' 

Therefore, u> can be determined from the 
measurement of the B° — B oscillation am- 
plitude: 

A m ix = P ° F : P J? F = (l-2w)cos(Am d At). 

FOF + FSF 

(4) 

The flavor of one of the two neutral B mesons 
is identified using B° — ► D*~ £ + v where D*~ 
decays to D tt~ followed by D decays to 
either K + tt~ , K + ir~n a , or K + n + Tr~Tr~ , or 
using B° — > D~£ + u where D~ decays to 
K + tt~it~ . We then identify the flavor of the 
accompanying B meson by applying the tag- 
ging methods described above. The vertex 
position of D*£v is determined by requiring 
the £ track and the reconstructed D momen- 
tum vector form a common vertex. The ver- 



Table 3. Summary of tagged events 



fcp 


N ev 


J 1 'tp K s {ir + n~ ) 


40 


J/^K s (n°n°) 


4 


ifj(2S)K s (TT+n~) 


2 


ifj(2S)(J/ipn+Tr-)K s (Tr + 7r-) 


3 


Xci{lJ/i>)K s (-K + K~) 


3 




4 




42 


Total 


98 



tex position of the tagging B meson is de- 
termined using the method described in the 
next section. The proper-time interval At of 
two neutral mesons is derived from the dis- 
tance between the two decay vertices. We 
obtain lj and Amj from a fit to the At dis- 
tributions of the OF and SF events with the 
expected functions, which include the effects 
of the At resolution and background. We de- 
termine Am d = 0.488 ± 0.026 ps -1 , which is 
in good agreement with the world average!. 
Figure [5] shows the measured A m i X distribu- 
tions together with the fitted functions. 

Table I summarizes the tagging efficiency 
and the wrong tag fractions. The total tag- 
ging efficiency is measured to be 0.52, and 
the total effective tagging efficiency (e e //), 
defined as a sum of e(l — 2ui) 2 over all tag- 
ging methods, is 0.22. The statistical error on 
sin2</>i is proportional to 1 /^/e e //- Table 
lists the number of tagged events for each 
fcp- We find a total of 98 tagged events, 
of which 14 events were tagged by a high mo- 
mentum e, 12 by a high momentum fi, 48 by 
K ± , 3 by a medium momentum e, 3 by a 
medium momentum /i, and 18 by a soft tt. 

5 Proper-time Interval 
Reconstruction 

The fcp vertex is determined with the two 
lepton tracks from J/ip. We require that at 
least one of two tracks must have SVD hits 
and that the vertex point is consistent with 
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Figure 6. The average shape of the event-by-event 
resolution function obtained by summing over 300 
B ± -+ J/il)K ± events. 

the IP profile. The ftag vertex is formed from 
the remaining tracks in the event. In order 
to reduce the bias due to long-lived particles, 
tracks are excluded from the fit: if when com- 
bined with any other charged track they form 
an invariant mass consistent with a Ks', if 
the track has a large tracking error in z di- 
rection (a z > 0.5 mm); or if the minimum 
distance between the track and the recon- 
structed fcp vertex is too large, Sz > 1.8 mm 
or 6r > 0.5 mm (in the r<j) plane). In addi- 
tion, the track with the worst x 2 is removed 
from the fit if the reduced x 2 °f the vertex 
fit is more than 20. This procedure is iter- 
ated until the reduced x 2 is below 20. The 
expected vertex resolutions are ~ 40 /im and 
~ 85 (im for the fcp and ftag vertices, re- 
spectively. 

The most probable At is estimated as 
Az//3jc. The resolution of At, R(At), is 
parameterized as a sum of two Gaussians, a 
main Gaussian arising from the intrinsic SVD 
resolution and the ftag vertex smearing due 
to the finite lifetime of secondary charmed 
mesons, and a tail Gaussian due to a few 
poorly measured tracks: 



The mean values (fx, fitail) an d widths (a, 
o'taii) of the two Gaussians are calculated 
event-by-event from the fcp and f tag vertex 
errors, taking into account the error due to 
the approximation of At ~ Az/(3"/c. The 
Gaussian parameters and f ta u(= 1 - fmain) 
are determined from the full MC simulation 
studies and a multi-parameter fit to B — > 
D*£is data. In addition, by measuring the 
lifetime of the D° — > K~tt + decays using 
only z coordinate information, we studied 
the intrinsic z vertex resolution. In order to 
show the average shape of the R(At), Fig. || 
was drawn by summing event-by-event R(At) 
functions over 300 B± — ► J/4>K ± (real) 
events. The width of R(At) is dominated by 
the main Gaussian (f ma in = 0.96 ±0.04); we 
find < a >~ 1.11 ps, < Gtaii >~ 2.24 ps, 
< (i >= —0.19 ps, where <> indicates the 
average over all events. (The /Mail was fixed 
to —1.25 ps based on the simulation stud- 
ies.) The non-zero negative mean values of 
the Gaussians reflect the bias in ftag vertex 
position due to secondary charmed mesons. 

Based on the above described vertex 
reconstruction and resolution function, we 
measured lifetimes of neutral and charged B 
mesonsa. Figure [?] shows At distributions 
of some of measured decay modes with the 
results of the lifetime fit. Table |4| summa- 
rizes the results. The obtained values for dif- 
ferent decay modes are consistent with each 
other and are in good agreement with the 
world averages!, t b o = 1.548 ± 0.032 ps, and 
t b ± = 1.653 ± 0.028 ps. This verifies the va- 
lidity of our Az measurement and R(At). 



6 Extraction of sin20i 

An unbinned maximum likelihood method is 
used to extract the best value for sin 2</>i. The 
probability density function expected for the 
signal distribution with a CP eigenvalue of 



Table 4. Summary of B meson lifetime measure- 
ments 



Decay mode 


lifetime (ps) 


B" -> D*+£~v 

B" -» D* + 7T- 
S -> D+7T- 

s° - J/^7T° 


1.50 ± omto 6 i 

-, rc+0.18+0.10 

J- 00 — 0.17-0.07 

1.4lig;^ ± 0.07 

i rc+0.22+0.09 
J- 00 — 0.19-0.15 


Combined 


1.50 ±0.05 ±0.07 


B U -> J/ipKs 
5° -» J/1,K L 


i r,+0.28+0.11 

i - d ^-0.24-0.19 
i qo+0.36 


B~ -> D*°£-V 
B~ -> D°tt- 
B~ -> J/tpK- 


1.54±0.10l l (J;o7 
1.73 ±0.10 ±0.09 

1 07+0. 13+0.07 
1 -° ' -0.12-0.14 


Combined 


1.70 ±0.061^0 



rjf is given by: 

Sig(At,r]f,q) = ^exp(-|At|/r B o) 
x{l — q(l — 2w)r)f sin 20i sin(AmdAi)} , 

(6) 

where q = 1(-1) if f tag = B°(B°) and u 
depends on the method of the flavor tag- 
ging as given m Table |. The values 
of t b o and Amid are fixed to the world 
averages!, 1.548 ± 0.032 ps and 0.472 ± 
0.017 ps -1 , respectively. By investigating 
events in background-dominated regions (the 
side bands in the AE vs M oeam scatterplot), 
we find that the probability density func- 
tion for background events for all fc p events 
(except for B° — > J/ipKz events) is consis- 
tent with Bkg(At) — exp(— \At\/T bkg ), 
where T k g = 0.73 ± 0.12 ps. The likelihood 
of an event, i, is calculated as: 

Pi = Vsig f*™ Si 9( S : Vf: ?)-R(At - s)d S 

+ (1 - Pstg ) /+~ Bkg(s)R(At - s)ds, 

where p S i g is the probability for the event 
being a signal, and R(At) is the resolu- 
tion function described in the previous sec- 
tion. The log-likelihood — In pi is calcu- 

i 

lated by summing over all signal events. The 
most probable sin 20i value is found by scan- 
ning over sin 2^i values to minimize the log- 
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Table 5. Results of CP fit to control data. 
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Figure 7. At distributions and the results of the life- 
time fit for (a) B° -> J/ip~K*°, (b) B~ -> J/i>K~, 
(c) B° -» J/V>A"s and (d) B° -4 J/ipK L . The lower 
solid curve represents the background distribution. 



Decay mode 


Apparent sin20i 


B u -► J/i/j(K+TT-)* a 
B~ -> J/V»if~ 

B° -> D*~l+v 


-0 094+ 492 
u.u»^_ 458 

i-U.^lO_ 238 

-0.096 ±0.174 
+0.09 ±0.18 



Table 6. Results of CP fit to tagged fcp events. 



Decay mode 


sin 20i 


J/iJjKs(7r + TT~) only 
All CP - 1 modes 
All CP + 1 modes 


+o-8ii8:SS 


All combined 


+0.45l^;S 



likelihood function. 

To test for possible bias in the analy- 
sis, we apply the same analysis program in- 
cluding tagging, vertexing and log-likclihood 
fitting to control data samples with null in- 
trinsic asymmetry: B>° — > J/if)K*°(K*° — > 
K+7T-); B~ -> J/ipK-; B~ -> D°tt-; and 
£?° — ► D*~£ + i> decays. Figure || shows the 
At distributions of B° -► J/tpK*°(K*° -> 
iT+TT-), B" -> J/ipK-, and B~ -> L>°vr- 
events. The results of the fit for apparent 
CP asymmetry are given in Table [j]; they are 
all consistent with null asymmetry. The 
results of the fit to the tagged fcp events 
are summarized in Table [| To display the 
fitted results, the dN/dAt distribution for 
q = +1 events and dN/d(—At) distribution 
for q = — 1 events are added: 

dN/dAt\ q=+1 +dN/d(-At)\ q =- 1 
oc exp(— \At\/TB») 

x{l — (1 — 2w)r)f sin 2^1 sm(AmdAt)} . 

(8) 

Figures ^ (a) and (b) show the results for 
only fcp — J/iPKs(tt + tt~) events and for all 
CP — 1 events combined. Figure || (c) shows 
the result for CP + 1 events, i.e. J/iJjKl 
and J/ipn°. In fitting to fcp = J/^K^ 
events, the background due to J /^K*{Klt: ) 
+ non-resonant J/ipKLir , which amounts to 
~ 17% of the total background, is taken to 
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Figure 9. dN/dAt\ q=+1 +dN/d(-At)\ q= -i distribu- 
tions and the results of the CP fit for (a) B° — > 
J/ipKg(ir+ir~) only, (b) all CP — I fcp modes com- 
bined, and (c) all CP + 1 fcp modes combined. A 
lower solid line in each figure indicates background 
contribution. 



Figure 8. At distributions and the results of the CP 
fit for (a) B° -> J/4>K*°(K*° -> K+n~), (b) B~ -> 
J/ipK-, and (c) B~ -> L> tt- events. 
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Figure 10. The log-likelihood as a function of sin 2<f>\ 
(a) and 
a sum of dN/d(+At)\ q=+1 +dJV/d(-At)| 9 =-i for 
•qj = -1 and dN/d(-At)\ q =+i +dN/d(+At)\ q =-i 
for tjf = +1 (b) for combined CP - 1 and CP + 1 
events. 



be a mixture of CP — 1 (73%) and CP + 
1 (27%) states, based jon the results of a 
B -> J/tpK S Tr° analysis!. A fit to 52 events 
in the J^jKl sideband (i.e. 1.0 < p* B < 
2.0 GeV/c region), where the non-CP J/ipX 
events dominate, gives the result sin 2<pi — 
+0.02lQ4g, consistent with null asymmetry. 
Finally we perform the simultaneous fit to 
CP — 1 and CP + 1 events to extract the 
best sin 201 value. Figure |l^ (a) shows the 
log- likelihood as a function of sin20i for a 
total of 98 CP - 1 and CP + 1 combined 
events together with the results for CP — 1 
and CP + 1 separately. We find sin 2^1 = 
+0.45±^. To display the results of 
the fit, dN/d(+At)\ q = +1 + dN/d(-At)\ q =-i 
for rj f = -1 and dN/d(-At)\ q = +1 + 
dN/d(+At)\ q= ^ 1 for Vf = + 1 ar e added 
so that the distribution becomes approxi- 
mately proportional to exp(— |At|/rBo){l + 
(1 — 2u>) sin 20i sin(ATO,iAt)}, as shown in 
Fig. 0(b). 

We generated 1000 toy Monte Carlo ex- 
periments with the same number of tagged 
CP events, having the same composition 
of the tags and the same resolutions as in 
the CP data sample, for an input value of 
sin 20i = 0.45. Figures |ll| (a) and (b) show 
the distributions of the central sin20i value 
and the statistical errors, + side and — side 
separately. We found that the probability 
of obtaining a value of the statistical error 
greater than the observed value was ~ 5%. 

Table lists systematic errors. The 
largest error is due to uncertainty in the 
wrong tag fraction determination. This was 
studied by varying to individually for each 
tagging method. The effect due to uncer- 
tainty in At resolutions for both signal and 
background was studied by varying parame- 
ters in R(At). Also included are the effects 
due to uncertainties in estimate of the back- 
ground fraction, in the world average tb and 
rrid values. An imperfect knowledge of the 
event- by- event IP profile could cause a sys- 
tematic error in sin 2<pi due to the vertex re- 
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Figure 11. Results of 1000 toy Monte Carlo exper- 
iments generated for sin2</>i = 0.45: distributions 
of the fitted sin20i (a) and statistical errors (b). 
In (b) distributions for errors on + side and — side 
are shown separately. Each experiment contains 98 
events with the same tag composition as the real 
data. 



Table 7. List of systematic errors of sin 2<f>i 



Source 




(7 — 


Wrong tag 


0.050 


-0.066 


R(At) 


0.026 


-0.025 


Background shape 


0.029 


-0.042 


Background fraction 


0.029 


-0.032 


t b o, Am d 


0.005 


-0.006 


IP profile 


0.004 


-0.000 


Total 


+0.07 


-0.09 




-0.004 -0.002 0.000 0.002 0.004 0.006 0.008 0.010 
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Figure 12. The region in the p — rj plane correspond- 
ing to this measurement ±1<t, sin20i = 0.45^p'^g. 
This was made by modifying Fig. 11.2 of ref. [5]. 



construction. The effect was studied by re- 
peating the entire fitting procedure by vary- 
ing the IP envelope by ±1ct in all 3 dimen- 
sions. The total systematic error of sin20i is 
found to be +0.07 - 0.09. 



7 Conclusion 

Based on a 6.2 fb _1 data sample col- 
lected at the T(45) resonance during the 
first year of the KEKB operation, we 
made a preliminary determination of sin 2<j)\ 
using 98 flavor-tagged events consisting 
of 40 J/ijKs(Tr+7r-), 4 J/^s(7r°7r°), 5 
^(2S)K s (7r+n-), 3 X ciK s (tt + tt-), 4 J/^tt 
and 42 J/i/jKl events. We found 

sin 20! = 0.45±°;g(stat)t°;°£(syst). 



Figure [12] shows the region in the p—rj plane^ 
corresponding to this measurement ±lcr, 
sin 20i = 0.45tg'45, together with the con- 
straints derived from other measurements!. 
While the current statistical uncertainty does 
not allow anything conclusive, this prelimi- 
nary result is consistent with the Standard 
Model prediction. 
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